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ABSTRACT 

The fracture problem of laminated plates which consist of ortho- 
tropic layers is considered. The orthotropic' material is assvimed to be 
of type II. Symmetrical cracks are located normal to the bimaterial 
interfaces. The external loads are applied away from the crack region. 
Three cases are considered; 

a) the case of internal cracks 

b) the case of broken laminates 

c) the case of a crack crossing the interface. 

A general formulation of the problem is given for plane strain and gen- 
eralized plane stress cases. The singular behavior of stresses at the 
crack tips and at the interfaces is studied. In each case the stress 
intensity factors are computed for various crack geometries. 

1. INTRODUCTION 

In recent years composites have attracted considerable attention 
largely because of their favorable crack propagation characteristics 
and of the flexibility they offer in designing a variety of structural 
components. In sheet structures experiments also show that the use of 
buffer strips with a relatively low stiffness may help to arrest a prop- 
agating crack. The composite materials are in general anisotropic. 
However, because of mathematical difficulties in most of the recent 
studies relating to composite laminates, the materials are assxamed to 
be isotropic [1-5]. An orthotropic strip containing a crack and bonded 
to two or thotropic half planes were considered in [6,?]. The problem of 
periodically arranged orthotropic strips containing cracks has recently 
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been investigated by Delale and Erdogan [s]. In this work it has been 
shown that regarding the roots of the characteristic eqviation orthotropic 
materials can be classified as of type I or of type II. In [s] the ma- 
terials of both layers are assumed to be of type I. The study of the 
elastic properties of orthotropic materials shows that, materials of 
type II are as practical as of type I. In the present study it is as- 
sumed that the composite plate consists of periodically arranged ortho- 
tropic strips of type II. Symmetrical cracks are located normal to the 
bimaterial interfaces and the external loads are applied away from the 
crack region. 

A general formulation of the problem is given for plane strain and 
generalized plane stress cases by the use of Fourier Integral Transform 
Technique. The resulting singular integral equations are solved numer- 
ically and the stress intensity factors are computed for various crack 
geometries. 

2. DISPLACEMENT AND STRESS FIELDS FOR STRIPS 


The composite plate shown in Figure 1 consists of two sets of 
periodically arranged strips of width 2hj^ and 2h2. Symmetrical cracks 
of length 2a and 2b are located normal to the biraaterial interfaces. 
The external load is applied away from the crack region. Using the 
superposition technique the solution of the actual traction-free crack 
problem may be obtained by superposing the solution of an uncracked 
homogeneous strip and the solution of a cracked strip where the self- 
equilibrating crack surface tractions are the only external forces 
(see Figure 2). For an ortho tropic solid the displacements satisfy 
the following field equations: 
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( 2 . 1 ) 


where 3j^(i=l/2f3) are elastic constants defined in Appendix A. The 
const-ants 3^ have different expressions for plane strain and generalized 
plane stress cases. In our analysis, we will consider the generalized 
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plane stress case. But the analysis is also valid for the plane strain 
case, simply by redefining the elastic constants. 

In [s], it has been shown that solutions satisfying equations ( 2 . 1 ) , 
yield the following characteristic equation: 


s“ + 3 ^s^ + 35=0 


The roots of (2.2) are: 


Sf = (0i + = / (-3^+3^) /2 

s^ = 03^ + iw^ = /(-34-3g)/2 


( 2 . 2 ) 


(2.3) 


s 
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-s 
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s and s are both real or complex conjugates. 

The material is defined as type I when both roots are real and as 
type II when they are complex conjugates. In [b] it has been assumed 
that the material of both strips is of type I. In this study, the 
orthotropic materials will be assumed to be of type II. For materials 
of type II the roots of equation (2.2) can be written as: 


S, = 0) + 10)^ 

i O 2 


S = 0) - 103 

2 o 2 


where o) is taken as positive. Furthermore noting that 
o 


u(x,y) = -u(-x,y) ; v(x,y) = -v(x,-y) 


and u and v vanish when y goes to infinity, for y> 0 , the displacement 
expressions can be obtained by adding the solutions defined by equa- 
tions (3.1) and (3.6) of [s]. 
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Thus, for materials of type II the displacements for each strip 
or layer are: 


-oo 

u(x,y) = — [a(c») cos (fjL)„ax) sinh(O) ax) 

nj 2 o 


+ c (a) sin (w^ax) cosh (co^ax) ]cosayda 


I 

^ [E(a)cos(o)20fy/»^,s) e 


+ G(a) sinCw^ay)//^) e 


]sinaxda 


v(x,y) = - 


{a (a) [g ' cos (co.ax) cosh (w ax) 

/ 2 O 

- 3 ^ sin(o)2CXx) sinh(w^ax) ] + c(a)[g!^ cos(a)20ix)* 

• cosh(O)^ax) + 3 ^ sin(o)20tx) sinh(o)^ax) ]}sinayda 

-00 

+ ^ {E(a)[-3' cos(a)2ay/'^)-3g sin(w 2 ay//^) ] 

■'o 

+ G(a)[3g cos (o) 2 ay) /v^) -3g sin(W 2 ay//^) ]} • 
-o)^ay//37 

• e cosaxda 


(2.4a,b) 


Differentiating (2.4a,b) and using the stress-strain relations, the 
stress expressions for generalized plane stress case are obtained as: 


Tr(l-V V ) 

0^(x,y) = j {2A(a) [-A 2 sin{ 0 ) 2 ax) 


2E 


• sinh(o) ax) + A, cos(a)_ax) cosh(u ax)] 
o 1 z o 

+ 2C(a) [A, sin(a).ax) sinh(o) ax) 

1 2 o 


+ A^cos (W-ax) coshCoj ax)] } acosayda 
2 2 o 
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+ Jo^E (a) [A2cos(W20ty /^) + A^sinCoj^ay//^) ] 

+ G(a) [-A^cos(a)2ay/-/3^) + A^sindo^oty/v^)] } • 
-w^ay/v^ 

• e acosaxda 


tt( 1 -v V ) 

— a (X.y, = 

« T •* 


{2A(a) [-A-sin(a)-Ox) • 

b 2 


• sinh(to ax) + A_cos (aj.ax) cosh(o) ax)] 
o 5 2 o 


+ 2C(a) [a sin{a)„ax) sinh(o) ax) 
b <6 o 


+ A^cos(a) ax)cosh(oj ax)] }acosayda 


^{E(a) [A^cos(o)2ay/>/g7) + AgSin (w^ay//^) ] 


+ G(a) [-AgCos( 0 ) 2 ay//^) + A^sinCw^^y/*^) 


-o)^ay//37 

e acosaxda 


•rr 


2g xy 
xy 


(x,y) = [ {2A(a)[A 


„cos(Ci) ax)sinh((D ax) 
» 2 o 


- A sin ( 0 ) ax) cosh (w ax)] + 2C(a) [A sin(a) ax) 

XU ^ O y 2 


cosh(0J ax) + A,^cos{(jJ^ax) sinh(O) ax) ]}asinayda 
o 10 2 o 

I {E{a)[-A^^sin (a)2ay/v^) + A^^cos (w^ay//^) ] 


+ G(a)[A^^cos(a)2ay/V^) + A^^sin (W2ay/v^) ]) 

V 

-w ay//^ 


• e 


asinoxda 


3^,3j and Aj^ are elastic constants given in Appendix A. A superscript * 
will be used for the elastic constants and the unknown functions when 
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these expressions are used for the second strip. 


3. FORMULATION OF THE PROBLEM 

The unknown functions A(a) , C(a) etc., which appear in the dis- 

placement and stress expressions can be determined by 
following boundary and continuity conditions; 

satisfying the 

Ui(hi,y) = U2(“h2/y) , 


Vi(hi,y) = V 2 (-h 2 ,y) , ( 0 <y<~) 

(3.1a,b) 



=■ ■'2xy<-'’2'y> ' • 

(3.2a,b) 

u^(0,y) = 0 , ^ ° ' (0<y<“>) . 

(3.3a,b) 

U2(0,y) = 0 , '^2xy^°'^^ " ° • 

(3.4a,b) 


^lxy<=‘l'‘» " ° - 


<hj^ , 


T (x„,0) - 0 , 

|x- i 

<h, . 

(3.5a,b) 

2xy 2' ■ 

' 2' 

2 


“ly<>‘l'<» = -Pl<V 

f 

lxj_|<a , 


Vi<Xi,0) = 0 , a< 

1x^1 

A 

b* 

H 

• 

(3.6a, b) 

OjyOCj.O) - 

t 

Ixjhb , 


^2^*2'°^ =0 , b< 

1*2* 

!<h2 . 

(3. 7a,b) 

conditions (3.3a,b) 

and 

( 3. 4a, b) are satisfied identically. 

Using 


(3.5a,b) we obtain 

^17 * ^12 * 

G(a) = - E(a) and G (a) = - ttf- E (a) , (3.8a,b) 

^11 ^11 
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Defining 


3Vi(Xi,0) 

9x, 


= ({)(x^) such that <f>(x^) = 0 for |x^|>a 


and 


3v2(x2,0) * * 

= (|) (x^) such that <{) (x^) = 0 for lx 2 Pb , 


the mixed boundazry conditions (3.6a,b) and (3.7a,b) reduce to the follow 
ing singular integral equations: 

-a (j,(t) 

dt + 

^o 


14 


<J>(t) 


{2A(a) [-AgSinCoj^ax^) 


sinh (o) ax, ) + A_cos (aj_ax, ) cosh(o) ax, ) ] 
o 1 5 2 i o ± 

+ 2C (a) [A^sin ((jo_ax. ) sinh(o) ax. ) 

5 2 1 o 1 


+ A_cos (o) ax, ) cosh{w ax, ) ]}ada 

6 2 1 Ox 


2E 


XV yx , , ^ 


and 


14 


b (t) 


^b “■^2 


dt + 


{2A* (a) [-A^sin(W 2 CtX 2 ) 


* * * * 

sinh(w ax„) + A_cos (o)„ax )cosh((jD ax ) J 

^ ^ ^ M O JL 


* _ * * * 

+ 2C (a)LA_sin(w^ax„)sinh{a) ax ) 

S 2 2 0 2 

.+ A-cos(a)^ax„)cosh(Q) ax )]}ada 
6 2 2 0 2 

* * 

7r(i-v V. J 

P 2 (x 2 ) -b<X 2 <b (3.9a,b) 

2E* 

y ■ 

Using the continuity conditions (3.1a,b) and (3.2a, b), the unknown 
functions A (a) C (a) , A* (a) and C* (a) are found to be: 
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A (a) = ^ [R^f^ + + R 3 h^ + R^m^] 

o 

c(a) = i [Rjfj + R^gj -f Rjhj + R^m^] 

o 

A (a) == -^ [Rj^fa + R^g^ + R^h^ + ^4”'3l 
o 

C*(0) = i [R^f^ + R^g^ + R3h^ + R^m^] 


(3.10a,b,c,d) 


The expressions for A^(a), R^(a), g^(a), h^(a) , m^(a) , (i=l,4) are 

given in Appendix B. 

Substituting (3.10a, b,c,d) into (3.9a,b), we obtain: 
ra 4>(t) ra rb 


ra ra m 

-a 1 i-, i-h 


^12 (t)dt 


TTCl-V V ) 
xy yx 

2E A, . 
y 14 


Pl(Xi) -a<x^<a 


fb <p (t) 

f 

— — dt + 

L ’=-*2 J 

k2i(x2,t)(j>(t)dt + 
-a 


(t)dt 


* * 

7r(l-v V ) 
xy yx 

y 14 


PjCXj) -b<Xj<b 


where 


-0) a(h,-t) 


f -u) wm, 

^11 ^^1'^^ ~ I 
.00 


da 


-w a(h -t) 
a)e da 


^21 ^^2 


"^22<^2 


»t) = I 

ft) = j 


° -w a(h -t) 

k 3 (x 2 ,t,a)e ° da 


’ -WQa(h2-t) 

k^(x2,t,a)e dot 


(3.11a,b) 


The expressions for (i=l,4) are given in Appendix B. 
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The integrands of the kernels (i,j=l,2) are bounded for all 
values of a, but are singular of order 1/a when a=0. In [s], it has 
been shown that using the single-valuedness conditions /^<}) (t) dt==0 cind 
(t)dt=0 the kernels k^^ can be integrated numerically. 

4. CASE OF BROKEN LAMINATES 

When one of the cracks touches the interface (i.e. , a=h^^ or 
we obtain the case of broken laminates. In this case, for example, 
for a=hj^ the kernel kj^j^ becomes unbounded as (Xj^,t)->ih]^ simultaneously. 
Therefore the kernel k^^^ must be considered in two parts: 


’'ll'*!'*’ ° ’'lls**!''"’ ’'llf**!''"’ 


where is the singular part and is the bounded part. 

After some lengthy algebra, the singular part of k^^^^ is found to 

be: 

k (x ,t) = F (x ,t) + F„(“X-,t) 

11s 1 s X S 1 


where 


F^(x, ,t) = 
s 1 


U)2(t-Xi) 


\q 1 63 Loj^ a-X^) + [w^(2h^-t-Xj^) Y‘ 


+ \ 


a32(2h^-x^-t) 


64 Lw2(2hj^-Xj^-t)]2 + [03^ {2h^-t-x^) J2 


0) (2h, -t-x, ) 

O 1 1 


^65 Lw2(t-x^) J2 + l.a)^(2h^-t-Xj^) J2 
W (2h -t-x.) 

O X X ■ .1 

66 [W2(2hj^-Xj^-t) J2 + Lw^(2hj^-t-Xj^) J 


+ A. 


(4.1) 


The elastic constants are defined in Appendix A. The governing sing- 
ular integral equations then become: 


ihi 
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• (j)(t)dt + 


L ’'12<*1' 


* TT(1-V V ) 

t)4, (Odt = - - 3E ~ y ^ Pi<Xi) -h^<Xi<hi 

y 14 


and 


/•b (J) (t) |-hi 


dt + I k2^(x2,t)ct)(t)dt + 1 k^^(x^,t)<l) (t)dt 


22 ' 2 ‘ 


TT(1-V* V* ) 

( x^ ) -b<x^ <b 


2E A 
y 14 


2 2 ' 


(4.2a,b) 


In equation (4.2b) since the only singular term is •— — , the singularity 

U.-X2 

power is 1/2 at the tip of the crack in the second layer. But in equa- 
tion (4.2a) we have further singular contribution from which may 

result in a singularity power diffei'ent than 0.5. To find this singular- 
ity power y, let us first write equation (4.2a) in the form; 


*-h, 1 


(4.3) 


where ^ bounded function for all values of x^. 

F ( t) 

Assuming (J)(t) = Y following the procedure described in 

[s], the following characteristic equation is found from which the sing- 
ularity power y can be determined: 


H\Y) = COSTTf + 


^67 X *"X ^2 

— t: — 5T"^ — cos(2Ytan — ) 

hs '“^“2”‘58 “o 




<“^“ 2>*-58 


-1 ^^2 

sln(2Ytan [%^63+V65^ ' ° 


(4.4) 


For practical orthotropic materials equation (4.4) has only one root 
between 0 and 1. 

5. CASE OF A CRACK CROSSING THE INTERFACE 

To formulate this problem we shall first consider the crack config- 
uration shown in Figure 3. Noting that. 
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* 

j (t)dt = j (t)dt (i=l,2) 


rb (}) (t) r® 1 1 * 

lb t-x, ~ Jo t+x„^ ^ 


J-b ^ -^2 2 ‘' *2 

the governing singular integral equations for this case can be written 


ra (j)(t) ra. rd ^ 

■ ^Vx k^^(x^,t)c()(t)dt + 

-a 1 ia ° 

Tr(l-v V ) 



y 14 




'21 '"2 ' 


•rr(i-\) V* ) 


■D ^ TUl-V V ; 

[k (x ,t) -k-„(x ,-t)]cf) (t)dt = ■ c<x <d 

J, 2 /- 2 22 2 ow A 2 


2 E A, . 
y 14 


If a=h, and d=h_ all the kernels k..(i,j=l, 2 ) in ( 5 . 1 ) become unbounded 
1 2 1 ] 

for (Xj^,t)=+hj^ and (X2#t)=h2 simultaneously. Therefore, the singular 
part of kernels k^^ must be separated. The kernels can be written as: 


k. . = k. . + k. . 

13 13 s i3f 


where k. . is the singular and k. is the bounded part of k. .. The 
1.3s ^ 13 f ^ 13 

expressions for the singular part of each kernel are: 


= E (x ,t) + F (-X ,t) 
11s 1 s 1 s 1 


k (x,,t) = G^(x ,t) + G^(-x ,t) 
12s 1 si s i 
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where 


Fs(Xi,t) is defined by equation (4.1) 


and 




58 


0)2 (hj^-x^) - 0)2 *h2-t) 


[ +[u)*(h2-t)+0)^(h^-x^)] 




+ X 


0)2 (h^-x^) + 0)2 (h2~t) 


82 


+ 


+ X, 


] +[o)*(h2-t)+U)^(h^-x^) ] 

0) (h,. -t) +0) (h -X ) 
0 2 oil 


[Wq ■*■“0 ^ ^^ 1 "^ 

* 2 
c^)-0)2(h2-t)] 

0) (h,-t) +0) (h -X ) 
02 oil 

—2 

* 2 


H^(x2,t) = 


68 


^2 ^^2”*2^ ”^2 


[o)2(h2-X2) -0)2(h^-t) ] +[o)^(hj^-t)+o)^(h2“X2) ] 


+ X 


^2 ^^2"^2^ ^^2 


74 r. * 




+ X 


0) (h. -t) +0) (h_-x ) 
o 1 o 2 2 


75 


[a)^(hj^-t) +u" (hj-Xj) ]%[u* (hj-x^) -Uj (hj^-t) ] 


+ X 


0 ) (h. -t) + 0 )^ (h^-x_) 

O 1 O 2 2 
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2 ^ * 


[o)^(h^-t)+0)Q(h2-X2) ] +[o)2(h2-X2)+0)2(h^-t) ] 



(I) (t-x ) 

I (X .t) r— s ; 

= = ^8 [ [a.2(t-X2)]%[a)^(2h2-t-Xj)]' 


+ X 




90 r* -Zf* .2 

[w2(2h2“X2-t) ] +[o)^(2h2-t-X2) ] 


+ X 


(U^(2h2-t-X2) 


91 p*- T^r* n2 

[a)2(t-X2)] +[o)^(2h2-t-X2)] 


+ X 


w (2h--t-x») 
O 2 2 


[a.* { 2h2-t-Xj) ] %[10* ( 2h2-X2-t) f 


(5. 2a,b,c) 


The singular integral equations then become: 

r 1 T 

— h- 1 — hn 


rh ^ fh- 

+ j ^123^^1'^^'^ (t)dt + [k^2f 

’’ c ^ ^ 


1T(1-V V ) 


2E A 


P-, (X, ) -h, <x, <h 


y 14 




1 1 1 


and 


f^2 pi In* r^l 

Jc t X2 t:+X2 


^213^^2 


f^l 

+ 1 k2j^^(x2,t)4)(t)dt + 

^hi 


h2 


k22s(x2ft)(|) (t)dt 


* * 

f^2 * 7T(1-V V ) 

+ [k (X ,t)-k (X ,-t)](t, (t)dt = - s2ps!_p (X ) 

^ 'Vl4 


c<x <h, 
2 £ 


(5.3a,b) 
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To find the singularity power at the interface, assume: 




(h^-t^) ^ 


(t-c) ^(h 2 ~t) ^ 


Considering the singular parts only, equations (5.3a,b) can be written 


as! 


1 r 1 T 1 * 

i’ J + k^j^(x^,t)]4.(t)dt + - , (t)dt 

-h, 1 ^ 


- Qj(Xj) 


and 


1 f^l 

IT J 
-h. 




h2 




- QjCXj) 0<Xj<h2 


(5.4a,b) 


where ^2 ^^2^ bounded for all and ^ 2 ' 

Following the procedure described in [s], one can obtain the eharac 
teristic equations as: 


cotir6 = 0 


and 


1 

A (3) = [cosFg + -™ + Agrees (Xgg3) + AggSin(Agg3)J 
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^93 1 

[cos7t3 + 


^68 
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(5.5a,b) 


Equation (5.5a) gives the expected 1/2 singularity power at the crack 
tip. But equation (5i5b) is rather complicated and its handling needs 
care. In [s] it has been shown that for some orthotropic materials of 
type I there is no singularity at the interface. This is also a possi- 
bility for materials of type II. For practical materials equation 
(5.5b) gives a root between 0 and 1. In deriving the characteristic 
equation (5.5b) one may note that there is a relation between F(h^) and 
F*(h^). 


a^F(h^) + a2F (h^) = 0 


(5.6) 


where 


a, = - 


(2hj^) 


^ [cos7i 3 + + Xg^cos(X^^3) + Xo^sin(X^^3)] 


3 


68 


96" 


‘9 5' 


and 


a„ = - 


*2 *2 

<“o *“2 > 


*2 *2 

fW +0)_ 
o 2 

2 2 

03 +W 
O 2 


96" 




This is a condition to be used in obtaining the solution. 

6. THE SOLUTION AND THE RESULTS 
6.1 Case of Internal Cracks 

In this case the governing equations are (3.11a, b) with the single- 
valuedness conditions 


c{) (t) dt = 0 , 


7b 


4) (t)dt = 0 


-a 


-b 
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The singular integral equations can be solved in a straight-forward 
manner by using the Gauss-Chebyshev integration technique described 
in [ 9 ]. In fracture problems one is generally interested in the stress 
intensity factors which can be expressed in terms of the density func- 
tions F(t) and F*(t). The stress intensity factors are defined as 
follows : 


For a<h, ; k = lim »/2(x -a) O’, (x,,0) 

1 a x.->a 1 ly 1 

(6.1a,b) 

and for b<h_: k, = lira /2 (x^-b) 0_ (x-,0) 

2 b 2 2y 2 

By making use of (6.1a,b) and equations (3.11a,b) after some algebraic 
manipulations we obtain; 


k 

a 


-2 


A, 

14 y 

(1-V V ) 
xy yx 


F (1) 
o 


k 


b 


-2 


* 4r , 

A, .E 
14 y 

★ * 

(1-V V ) 
xy yx 


F*(l) 

o 


where the index "o" denotes the normalized quantities. The loads p^, 
p^ are not' independent. Assuming that there is no constraint in x 
direction, p^ and p^ satisfy the following relation: 



In the computation the following material combination is used: 


Material of the first strip: 

E = 3.1x10^ psi (21.3745x10^ N/m^) 

E = 9.7x10^ psi (66.8815 x10^ N/m^) 

y 

G = 2.6x10^ psi (17.927x10® N/m^) 
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0.200 


V 

xy 

Material of the second strip: 

= 2.5x10^ psi (17.2375 x 10^ N/m^) 

Ey = 2.5x10® psi (17.2375x10^ N/m^) 

= 1x10^ psi (6.895x10^ N/m^) 


V = 0.760 
xy 


By letting a, b, or go to proper limits or choosing the materials 
close to isotropic, one can recover all the special cases done in [5], 
[6] and [lO]. 

Figures 4 and 5 show some of the computed results. In Figure 4 
the stress intensity factor k is plotted versus h_/h, for the case 
b=0 (i.e., there is no crack in the second strip). For h2=0 one re- 
covers the results for coll inear cracks imbedded in a homogeneous me- 
dixim. When h -x», k reaches an asymptotic value which can be found in 

Si 

[6]. Figure 5 shows the stress intensity factor kj^ for the case when 
a=0 (i.e., there is no crack in the first layer). 

6.2 Case of Broken Laminates 


In this case the solution may be obtained by considering equations 
(4.2a,b) and the single-valuedness conditions 


rhi 

(j) (t) dt = 0 




(j>* (t)dt = 0 


The singularity power y at the bimaterial interface is obtained from 
the characteristic equation (4.4). The singular integral equations 
are solved by using the numerical method described in [9]. For this 
case the Gauss-Jacobi integration method is used and the fundamental 
functions are the weights of Jacobi Polynomials (for more details see 
[s]) . The stress intensity factor is defined as: 
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i 


% . 


For a = h : k = Ixm 2^(x +h-)^a (x ,0) 

-■■a ^ z ^ 

2 "2 

Using definition (6.2) and equation (4.2b) k is found to be; 

& 


k = 
a 


* * 

2E A- . 

y 14 


V 


1c * 0 0 

(1-v V ) (co'^-:ujM 
xy yx o 2 


^2+0)^ 
^ o 


*2 *2 

(W- +W 
2 O 


Y 




( 6 . 2 ) 


-1 ^2 -1^2- 
• I cos[Y(tan ;;;^ + tan ^ 


-1 ^^2 -1 ^2 

+ sin[Y(tan — + tan — sr) ] (“^_^7-i“^9^7c) 

0) 0J„ o Z /D 


, V 00 

tO^ “1 2 T 

+ cos[y(-tan jj- + tan ^) ] (-W2^74-<^o^7e) 

O 


CO 

+ sin[y(-tan‘^ ^ + tan"^ •— -) ] (“<^o>^ 74 +f^ 2 ^ 75 > 


Figure 6 shows the variation of k when the first laminate is broken. 
The variation of k^^ when the second J.aminate is broken is shown in 
Figure 1. 

6.3 Case of a Crack Crossing the Interface 

For this case the equations to be solved are equations (5.3a,b), 
the single-valuedness condition /||^(})(t)dt=0 and relation (5.6) . The 
singularity power 3 at the interface is obtained by solving equation 
(5.5b) n\3merically. Using the Newton-Raphson method for the material 
combination considered it is found that 

3 = 0.0852. 

For the solution of the singular integral equations again the Gauss- 
Jacobi Integration Technique described in [o] is used. In this case 
the fundamental functions are the weights of Jacobi polynomials. The 
stress intensity factors can be defined as follows; 
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(6.3a) 


k = lifi (c-x ) a (x ,0) 
t> X2"^c ^ 2 


and at the bimaterial interfaces 


g 

k = lim Y' CT, (h, ,y) 

XX lx 1'^ 

g 

k = lim V T, (h. ,y) 

Y-V0+ " ^ 


(6. 3a,b) 


By making use of definitions (6.3a,b,c) and after lengthy algebra we obtain 


_ . * _* 
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2 simr3/2(ca^ ^\i 3 


-1^*2 

cos($tan ~-r ) + (-w_X, ,_+a)_X, - .) 




o 117 2 116' 


* ) 

-1 0)2 1 
sin(3tan — r ) ] V 


and 




k = 2G <, 
xy xy 


2^cos'ir3/2 


^25 ^2 

— ^ sin(gtan — ) 

2 U)o 




-1 ^2 
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7 ' 


-1 ^^2 1 

+ (X^^ga)2+X^lga)o)cos(3tan — ) 

o 


(>^) - c ) V ( l ) 

J 2 0 

3+*5 * 2 * 2 

2- \osiT3/2(a3o +W2 


^^“^12l“o'''\20‘^2^ 


-1 ^2 * * 
sin(3tan -r) + (X (0-+\._.a)^) • 

12X 2 X20 O 

* 

-1 ^ 

cos(3tan -|-) ] 


Figure 8 shows the variation of with c/h 2 , for different values of 

(hi/h 2 ) ratio, increases as (h^/h^) increases. Figures 9 and 10 

show the variation of k and k with respect to c/h«. 

XX 2 


ACKNOWLEDGEMENTS 


The author gratefully acknowledges the timely advice and guidance 
of Professor Pazil Erdogan. The support of Professors Erdogan and 
Arin is greatly appreciated. The work was supported by the Materials 
Division, NASA-Langley, under Grants NGR-39-007-011, NSG-1178 and by 
NSF under Grant ENG-73-045053 AOl. 


20 


APPENDIX A 


Definitions of the material constants: 
the material in the second strip. 


A superscript * will be used for 
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APPENDIX B 


Expressions of the functions used in equations (3.10a,b,c,d) 
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itZ #2 
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+ g*(a) (Aj^ 2 ^(a)h(a) + - X^^gia)h(a) + Xj^25(a)m(a) ) 

+ (-X,^h*(a) + X,^m*(a)) (g^(a) + f^(a))] 

1 / 16 

h,(a) - A, Jf*(a) (-X.g(a)h(a) + X f(a)m(a) 

3 10 o o 

- X^g(a)m(a) - X^f(a)h(a)) + g*(a) (X^f (a)h(a) 

+ X^f(ei)m(a) - X^g(a)h(a) + X^g(a)m(a)) 

5 -5 4 

+ {f^(a) + g^(a) ) (Xgm*(a) - Xgh*(a))] 
m^(a) = X239*(a) (h^(a) + m^(a)) 

+ h* (a) (X^^f (a)h(a) + - X^gg(a)h(a) 

+ X24g(a)m(a)) + m*(a) (X2gf(a)h(a) + X^^f{ct)m(a) 

- X^^g(a)h(a) + X 2 gg(a)m(a) ) 

f^(a) = (h^(a) + m^(a)) (X 2 gf*(ot) + X^^g (ot) ) 

+ h* ioL) (-X^^f{a)h{a) - X^^f{(x)m{a) 

+ X^^g(a)hia) - X^^g(,0L)m{a)) + m*(a) (X^^f (ot)h(a) 

+ X2Qf(a)m(a) - X2Qg(a)h(a) + X^^g(a)mm) 

g^(a) = [f*(a) (X^2 ^ + X^^f{a)m[CL) 

- X^^g{a)h(a) + X^^g{a)m{a)) 

+ g*{ot) (Xj^^g(ot)h(a) - X^^f(a)m(a) 

+ Xj^gg(a)m(a) + X^^^f (a)h(a) ) 
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h^(a) = -A^Q[f*(a) (A^f (a)h(a) + A5f(a)m(a) 

-A^g(a)h(a) + A^g(a)m(ay) + g (a) (A^g(a)h(a) 

-Agf(a)m(a) + A^g(a)m(a) + A^f(a)h(ot)) 

+ (f2(a) + g^(a)) (Agh*(a) + Agin* (a))] 

m^(a) = -Aggf*(a) (h^(a) + m^(a)) + h*(a)( 

Ag7g(a)h(a) - A^7f(a)m(a) - A^gg(a)m(a) 

- A^gf(a)h(a)) + ra*(a) (-A^^g(a)h(a) 

+ A2^f(a)m(a) + A^^g(a)m(a) + A^^f(a)h(a)) 

f(a) = cosCw^ah^) sinh(w^ah^) 

g(a) = sin (w_ ah- ) cosh(w ah.) 

2 1 o 1 

h(a) = cos (w„ah, ) cosh (oj ah.) 

2 1 o 1 

m(a) = sin(a)„ah_ ) sinh(oj ah^ ) 

2 1 o 1 

* * * 

f (a) = cos(o) 2 Cdh 2 ) sinh{a)^oh 2 ) 

lit "is 4 > 

g (a) = sin (oj^ah.) cosh (o) oh.) 

2 2 o 2 

* * * 

h (a) = cos (o)„ah„) cosh (o) oh„) 

2 2 o 2 

it it it 

m (a) = sinCw^oh^) sinhCoj^oh^) 



The expressions of k^(i=l,4) used in equations (3.11a,b) 

ki(Xi,t,a) = - — { 2 [-AgSin(W 20 ix^) sinh(to^ax^) 
14 o 

+ A_cos(co-ax, )cosh(u3 ax,)] * [f, (a)A, . • 

D Z 1 O J. 1 XCJ 

(A^^cosCw^aCh^-t) ] + sin[oJ2a(hj^-t) ]) 

+ g^(a) ^ (cos[o)2a(hj^-t) ] + A^gSin[u2a(h^-t) ]) 

+ hj^(a) A22(sin[w2a(h^-t) ]) 

+ m^(a) A^^ (— cos[w2a(h^-t) ] - sin[w2a(h^-t) ])] 

+ 2[A^sin(a)^ax, ) sinh(o) ax.) + A-cos(W-ax, ) • 

5 21 Ol 6 21 

cosh (co^ax^) ] • 

• [f^ (a) A^g(A^^co&[o) 2 a(h^-t) ] + sin[o32a(h^-t) ]) 

+ 92 ^“^ ’i’ (cos[o)2a(h^-t) ] + A^gSin[o)2a(h^-t) ]) 

+ h2(a) A22sin[w2a(h^-t}] + m2(a)A2^ ( 

— cos[w_a(h,-t) ] - sin[o) a(h -t) ])]} 

(Jj 2 1 2 X 

O 

■ A - 7A-(a T ^ ■ 

14 o 

sinh(03 ax, ) + A cos (w ax ) cosh(w ax, ) ] • 

O i. -> 4 X OX 

[f^{a) A^g(A^^cos[a)2 aCh^-t)] + sin[a32a(h2~t) ]) 

- g^(a) j (cosCu^aCh^-t)] + A^gSinCoj^a (h 2 “t)]) 
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-h^(a)A22A^sin[a)2“^^2”^^ ^ 

★ 

* ^2 * * 

+ ‘^^^[ aj ^ adi ^- t ) ] - sinCw ^ aCh ^- t ) ]) ] 

+ 2[A sinCw.ax, )sinh(w ax.) 
b 2 1 o 1 

lAf 4> 

+ A cos(o) ax )cosh(w ax ) ]• [f (a) A. . (A • 

6 2 1 o 1 2 16 17 

A * 

cos[o)2a(h2-t) ] + sin[o)^a(h2-t)]) 

- 92^“^ 4 (cos[o)2a{h2-t)] + A^^ sin[o)2a(h2-t) ]) 

* r * -1 

- h^Ca) A22A-j^sin[w2a(h2-t) J 

* 

* ^2 * 

- sin[o)2a(h2-t) ])]} 

■ A*J (a) ‘"['4 =“<“ 2 “ 2 > • 

±4 O 

* * * * _ 

sinh(w ox„) + A_ cos ((jO„ox„) cosh (w ox.) ] * 
o 2 5 2 2 o 2 

[f^Ca) A^g(A^^cos[o)2a(h^-t) ] + sin[o)2a(h^-t) ]) 

+ 93 (Oi) ■|•(cos[(J02a(h^-t) ] + A^gSin[o)2a(h^-t) ] ) 

+ h^(a)A22sin[w2a(h^“t)] 

+ m^(a) A2^(1^ cos[t02a(h^-t) ] - sin[oj2a(h^-t) ]) ] 

+ 2[Ag sin((j02OX2) sinh(a3Qax2) 

* * * 

+ A,cos(o)„ax„)cosh((o ax„)J • 

6 2 2 o 2 ■“ 
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• [f^ (Ot) A^g(A^^cos[t02a(h^-t) ] + sin[w2a(h^-t) ]) 

94(a) j (cosEw^aChj^-t) ] + Aj^gSintoj^adi^-t) ]) 

+ h^(a) A 22 sin[o) 2 a(h^-t) ] 

0^2 

+ ”'4^a)A24 <^:;^ cos[w2a(h^-t)] 

-• sinCw^ax'h^-t)])]} 

k^<Xj,t,a) = {2 [-A*sin(ai*ax2). 

it i( ic i( 

sinh(w ax„) + A_ cos (co^ax.) cosh (co ax„)l • 

O Z 5 22 0 2 

•it "it ’h 

• Ff2(a)A^g(A^^cos[a)2a(h2-t)] + sin[o)2a(h2 ' ^H) 

I r* T * r* 

- g^Xa) -4 (cosLa)20t (h2“t) J + A^gSin[o)2a (h2~t) J ) 

- h3(a)A*^X^sin[a)*o(h2 -t)] 

* 

* * 

+ m^(Ot) ^2^24 ^ 

^o 

- sin[o)*a(h2-t) ]) ] 

•it 4r 

+ 2[A sin(co ax„)sinh(co Ctoc„) 

5 2 2 o 2 

it it it 

+ A cos (0)^0ix„) cosh (0) Ox^) ] • 

6 2 2 O 2 

it it ' " it it 

[f^ (a) A^^CAj^^cosCm^aCh^-t) ] + sin[w2a(h2-t) ]) 

1 it it it ... 

"■ 94(a) J (cos[aJ2a(h2-t) 3 + A^gSin[w2a(h2-t) ]) 

- h4(a) A*2X^sin[w*a(h2 -t)] 

* 

* S * 

+ m^(a)X2A24 (-F cos[a32a(h2-t) ] 

^o 

- sin[a)*a(h2-t)])]} 
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Figure 4. The stress intensity factor for the crack in the first 
scrip. 







Figure 7 . The stress intensity factor kj^ when the second laminate is 
broken. 
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Figure 10. The stress intensity factor k^y for a crack crossing the 
interface. 


